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ABSTRACT 

We discuss in some detail the production of extremely high energy (EHE) neutrinos 
with energies above 10 18 eV. The most certain process for producing such neutrinos 
results from photopion production by EHE cosmic rays in the cosmic background 
photon field. However, using assumptions for the EHE cosmic ray source evolution 
which are consistent with results from the deep QSO survey in the radio and X-ray 
range, the resultant flux of neutrinos from this process is not strong enough for 
plausible detection. A measurable flux of EHE neutrinos may be present, however, if 
the highest energy cosmic rays which have recently been detected well beyond 10 20 eV 
are the result of the annihilation of topological defects which formed in the early 
universe. Neutrinos resulting from such decays reach energies of the grand unification 
(GUT) scale, and collisions of superhigh energy neutrinos with the cosmic background 
neutrinos initiate neutrino cascading which enhances the EHE neutrino flux at Earth. 
We have calculated the neutrino flux including this cascading effect for either massless 
or massive neutrinos and we find that these fluxes are conceivably detectable by air 
fluorescence detectors now in development. The neutrino- induced showers would be 
recognized by their starting deep in the atmosphere. We evaluate the feasibility of 
detecting EHE neutrinos this way using air fluorescence air shower detectors and derive 
the expected event rate. Other processes for producing deeply penetrating air showers 
constitute a negligible background. 

Subject headings: cosmic microwave backgrounds - cosmic rays - cosmic strings - 
early universe - elementary particles - instrumentation : detectors 

1. Introduction 

It is well known that there exist extremely high energy particles in the Universe with energies 
up to ~ 10 20 eV. How they are produced and how they propagate in the Universe have been 
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puzzles for a long time. Recent experiments suggest that these particles come from outside the 
Galaxy ( Bird et al. 1995 , Lawrence et al., 1991 , Yoshida et al., 1995| ). More detailed understanding 



of their origin will require detecting many more particles. It may also be necessary to observe the 
production process through other windows. Nucleons and nuclei lose their energies in less than 
lOOMpc through interactions with the cosmic microwave background ( preisen 1966 , Zatsepin 



and Kuzmin 1966| ). Furthermore, intergalactic magnetic fields can bend their trajectories, masking 



their sites of origin. Gamma rays are even more limited in pathlength due to electron pair 
production when they collide with radio photons. 

In this regard, neutrinos have uniquely advantageous characteristics: they can penetrate 
cosmological distances in the Universe and their trajectories are not deflected because they have no 
electric charge. They carry information about extremely high energy (EHE) production processes, 
even in the early Universe. It is therefore important to understand the possible processes for 
producing EHE neutrinos and to consider the possibilities for detecting the predicted flux resulting 
from several different models. 

What mechanisms might produce EHE neutrinos? It has been pointed out that they should 
be produced at least by the Greisen mechanism: the decay of photopions produced by EHE 



cosmic ray protons colliding with the cosmic thermal background photons (Greisen 1966, [Hil 



and Schramm 1985| , |Yoshida and Teshima 1993|) . The detection of the Greisen neutrinos would 
itself supply firm evidence that the EHE cosmic rays are coming from extragalactic space. The 
intensity of these neutrinos depends strongly on assumptions about the production of extragalactic 
EHE cosmic rays, and we need to make clear what kind of extragalactic sources might produce a 
detectable flux of Greisen neutrinos. 

Another possibility has been proposed which would result in more copious EHE neutrinos. If 
monopoles and/or cosmic strings (so called topological defects) were formed in symmetry-breaking 
phase transitions in the early universe, they may have produced EHE particles with energies up to 
the GUT scale (typically ~ 10 15 GeV) through their collapse or decay, with leptons and hadronic 
jets emitted from the supermassive "X" particles ( piiattacharjcc et al., 19*92] , Hill et al. 1987). 



This hypothesis is of interest here because these jets are expected to produce many more neutrinos 
than nucleons, and the EHE neutrino flux might be detectable. 

The detection of these EHE neutrinos is an interesting challenge. Currently there is no 
detector designed specifically for detecting EHE neutrinos. However, an EHE neutrino could 
induce an electron in the terrestrial atmosphere, thereby initiating an electromagnetic air shower. 
Several detectors with huge acceptance of ~ 10 4 km 2 sr for giant air showers are now being 
constructed and/or planned and we should not neglect the potential for EHE neutrino detection 
as a by-product of these detectors. The feasibility of the EHE neutrino detection by giant air 
shower detectors must be evaluated together with a study of the possible background event types. 



In this paper we analyze the production of EHE neutrinos and estimate their flux based 
on several different models. We also discuss prospects for EHE neutrino astrophysics which 
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would result from detecting such a flux. The paper is organized as follows: Section 2 reviews the 
Greisen neutrinos which are produced by photopion production of the EHE cosmic ray nucleons in 
intergalactic space. We present the possible flux for different assumptions about evolution of the 
EHE particle emitters and discuss the detection possibility. Some notes on neutrino emission from 
Active Galactic Nuclei (AGN) are given in section 3. In section 4 we focus on the emission of EHE 
neutrinos from collapse or annihilation of topological defects (TDs) such as monopoles, cosmic 
strings, etc. Because neutrino cascades with the cosmological relic neutrinos have a significant 
effect on the resultant flux in this scenario ( Yoshida 1994 ), we discuss the propagation of EHE 
neutrinos in detail for both massless and massive neutrinos in deriving the possible flux. We 
discuss in section 5 how to detect EHE neutrinos as deep air showers, using the air fluorescence 
detectors which are currently being constructed and planned. The event rate in these detectors 
for the TD scenario is also presented. The possibility of background events which could cause 
misidentification of neutrinos is considered in section 6. We summarize our conclusions in section 
7. 



2. The Greisen Neutrinos 

EHE neutrinos above 10 18 eV can be produced by photopion production between EHE cosmic 
ray nucleons and the 2.7 K background photons during propagation in intergalctic space ( [Hill and 



Schramm 1985| ). Because neutrinos travel much longer distance than the cosmic ray nucleons, 



the flux of these neutrinos depends heavily on assumptions about the evolution of the cosmic ray 



sources, including cosmic ray emission rates at high redshift epochs ( Yoshida and Teshima 1993 



Paper I). The source evolution function, r](t), is related to the cosmic ray spectrum as follows: 

J(E t0 ) = J'" dtMte)n(t e ) (f||y) 3 /^ dE te G(E te ,E t0 ,t e )f(E te ). (1) 

Here R(t e ) is the scale parameter of the Universe at time t e , n(t e ) is the number density of sources 
at time t e , Et e is the EHE cosmic ray energy at the emission time t e , Et is the energy after 
propagation, and f(E te ) is the energy spectrum at the source. G(E te , E to ,t e ) gives the energy 
distribution at the present epoch to for EHE cosmic rays which were input at time t e with energy 
Et e , resulting from energy loss due to interactions with 2.7 K photons and adiabatic loss due to the 
expansion of the Universe. This quantity plays the role of a "Green function" and is calculated by 
the transport equation for cosmic rays ( Hill and Schramm 1985| ) or Monte Carlo simulation for 



propagation in the 2.7 K photon field flYoshida and Teshima 1993| ). For astrophysical sources like 
radio galaxies or AGNs, we can consider the number of sources to be approximately conserved 
during the time scale we are interested in here. So Eq. (1) becomes 

J(E t0 ) = n f Z dt eV (t e ) f dE te G(E te ,E t0 ,t e )f(E te ), (2) 
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where no is the number density of sources at present. If we assume rj(t) to be proportional to 
(1 + z) m , where z is the value of redshift at time t, then Eq. (2) can be written as 



H, 



o Jo 



5 



J(E t0 )dE t0 = ^ I dz e (l + z e ) m -l I dE te G(E te ,E t0 ,z e )f(E te ). (3) 



Here Hq is the Hubble constant and we assume an Einstein-de Sitter Universe. In this expression 
m = corresponds to the case of no evolution. 

The neutrino flux becomes higher as the evolution parameter m and the "turn-on time" 
z max are increased. Figure [I] shows the spectrum of Greisen electron neutrinos under several 
assumptions for m and z max . These spectra are estimated by Monte Carlo simulation as in Paper 
I. The primary spectrum of cosmic ray nucleons is assumed to be ~ E~ 2 up to 10 22 eV though this 
assumption may be optimistic considering the proposed acceleration efficiency at radio galaxies 



( Biermann and Strittmatter 1987 ). However, even with these optimistic assumptions, the neutrino 
flux is found to be at most comparable with that of the EHE cosmic rays at around 10 19 eV. 
Although the neutrino flux would be a good probe of m and z max pertaining to evolution of the 
EHE cosmic ray production, detection of these fluxes with the current techniques is unlikely 
because of the low cross sections for neutrino interactions. 

Let us consider the case of extremely strong source evolution here. If the primary spectrum 
of EHE cosmic rays is very hard and emissions at high redshift epochs are responsible for the bulk 
of cosmic rays in intergalactic space, they would yield more neutrinos during their propagation. 
Figure || shows the expected neutrino spectrum for the case m = 7 and z max = 5 with a primary 
cosmic ray spectrum ~ E" 1 ' 5 calculated by the same Monte Carlo simulation as in Paper I. The 
neutrino flux exceeds the cosmic ray flux by a factor of ~ 140 at around 10 19 eV and may be 
detectable. However, we should note that these strong evolution assumptions would be inconsistent 
with the evolution of the luminosity functions of QSOs recently measured in the X-ray and radio 
range (Boyle et al., 1994, Dunlopc and Peacock 1990| ). The deep sky survey with the ROSAT 



X-ray satellite found that m = 3.25 and z max = 1.6 gives the best fit to the observed data (Boyle 



et al., 1994). Therefore if the powerful astrophysical sources like QSOs are the origin of the EHE 



cosmic rays as proposed by many authors (Rachen and Biermann 1993), it is hard to argue for a 
stronger evolution than this. 



3. The Diffuse Neutrinos from AGN 

Active Galactic Nuclei (AGN) are possible acceleration sites for producing EHE cosmic rays, 
and the accelerated proton energy loss due to pp and/or p7 interactions in the AGN accretion 
disk or with UV photons in the associated jets are the dominant mechanisms for neutrino (and 
photon) production. Many authors have studied these processes to predict the possible neutrino 
fluxes (e.g., (Szabo and Protheroe 1994 , Stecker et al. 1991 , Stecker et al. 1992] ) mainly focusing on 



TeV and PeV neutrinos. Whether AGN can produce EHE neutrinos with energies of 10 EeV or 
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greater depends on the proton maximum energy. Radio-quiet AGN, the most numerous class of 
AGN, are not likely to produce EHE cosmic rays while radio-loud AGN might be able to produce 
EHE cosmic rays and neutrinos. Mannheim modeled the production mechanism of EHE neutrinos 
in the jets of radio-loud AGN (Mannheim 1995) and found that the maximum neutrino energy 
reaches about 10 EeV. When the jet is pointing at the observer, particles accelerated in the 
moving jet plasma can gain an extra so-called Doppler factor ( Mannheim 1993 ). An uncertainty in 
this model arises from the fact that we don't know how many radio-loud AGNs contribute to the 
bulk of energetic neutrinos and photons and it would be difficult to give a reliable estimate of the 
normalization of the present-day neutrino flux. One possibility is to assume that the radio-loud 
AGNs are also responsible for the diffuse 7-ray background above 100 MeV (Mannheim 1995). 
The increasing number of EGRET sources and the detection of TeV photons from Mrk 421 and 
501 (punch 1992| , |Q uinn et al., 1995]) , support this hypothesis, though these photon emissions 
could be driven not by 7r° decay but by inverse Compton scattering of energetic electrons without 
producing any neutrinos. Under the assumption that the energy flux of 7-rays emitted from 
radio-loud AGN is comparable with the diffuse 7-ray background, one obtains the EHE neutrino 
flux as ~ 6 x 10 _33 m _2 sr _1 sec -1 eV~ l at 10 EeV ( Mannheim 1995) ) , which is higher than the 
cosmic ray flux by about a factor of 4. In this case the AGN neutrinos would dominate the Greisen 
neutrinos even in the EeV range. However, the predicted flux is still low and hard to detect, as 
will be described in the section 5, and detection of TeV-PeV neutrinos from AGN by underground 
neutrino telescopes is more feasible since the neutrino fluxes in this energy range are much higher 
and radio-quiet AGN can also produce TeV-PeV neutrinos. PeV neutrinos from AGN should be 
much easier to detect than EeV neutrinos due to the much higher expected fluxes. 



4. EHE Neutrino Emission from Topological Defects 



The possibility of producing EHE particles by the annihilation or collapse of topological defects 



(TDs) such as monopoles, cosmic strings, etc., has been proposed recently (e.g.,Bhattacharjee 
|et al, 1992| , |Hill et al. 1987) ). The maximum energy of the particles produced from TDs can 
reach the typical GUT energy scale. The energy trapped in the defects is released in the form of 
supermassive gauge bosons and Higgs bosons which are usually referred to as X-particles. The 
decay of the X-particles can give rise to quarks, gluons, leptons, etc., which materialize into EHE 
nucleons, photons, and neutrinos with energies up to the GUT scale. This mechanism could be the 
origin of the highest energy cosmic ray events well beyond 10 20 eV detected by the Fly's Eye ( [Bird 
et al. 1995| ) and the AG AS A ( Hayashida et al. 1994[) . Though we cannot exclude the possibility 
that these events could have originated at a nearby radio galaxy at an earlier time of strong 



activity if a strong magnetic field in the propagation space bent their trajectories ( Elbert and 



Sommers 1995| ), the TD hypothesis appears to be a possible option because these extremely high 
energies would result from the direct production of particles from the decay of the supermassive 
X-particles without any acceleration mechanism, and the lack of any identifiable astrophysical 
sources near the arrival directions of these events is not a problem in this scenario because TDs 
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are not necessarily associated with astrophysical sources ( Sigl et al. 1994| ). 



Among possible kinds of TDs, closed loops of cosmic strings have been studied most 
extensively as emitters of EHE cosmic rays. It has turned out that whether they would produce 
an observable cosmic ray flux depends on the abundance of loops in the lowest oscillation state 
which collapse within one oscillating period before losing their trapped energies by gravitational 
wave radiation. This is a free parameter, and no definitive theory gives a strict limit for EHE 
cosmic ray production by cosmic string loops. 

Recently the production of EHE cosmic rays through a process which involves formation 
of metastable monopole-antimonopole bound states (monopolonium) and their subsequent 
annihilation has been studied QBhattacharjee and Sigl 1995| ). In this scenario the measurable EHE 
cosmic ray flux requires a fractional abundance of monopolonium of more than ~ 1(T 7 of the 
total monopole abundance. Though this requirement is not excluded by the Parker bound, the 
monopolonium fraction is quite speculative. 

Therefore, the observation of EHE particles seems to be essential for testing whether these 
hypotheses are true or not. One prediction from the TD scenario is a high flux of EHE neutrinos, 
because the hadronic jets produced by the collapse of TDs would create a lot of pions which in 
turn emit neutrinos through their decays. The detection of the EHE neutrinos as well as the EHE 
gamma rays would support the hypothesis that the origin of EHE particles is related to GUTs. 

The energies of EHE neutrinos produced by the collapse of TDs range up to the GUT energy 
scale (~ 10 16 GeV) and their interactions with the cosmic thermal neutrinos (present temperature 



1.9 K) produce a significant neutrino cascade ( Yoshida 1994 ). The cascading increases the total 



EHE neutrino flux because a single neutrino from the X-decay can spawn multiple EHE neutrinos. 
Including this enhancement yields an EHE neutrino flux which is more likely to be detectable. In 
this chapter we derive the expected neutrino flux from TDs taking into account the enhancement 
factor due to neutrino cascades. Our treatment of the propagation of neutrinos in the cosmic 
massless neutrino background field is described elsewhere ( |Yoshida 1994 ), but we also present a 



brief description here because it is essential for the flux estimation. Then we further discuss the 
neutrino cascading effect for the the case of massive neutrino backgrounds. Finally we calculate 
the differential neutrino flux for several assumptions for the masses of "X" particles and neutrinos. 



4.1. Propagation of EHE neutrinos: Massless Neutrinos 

Because energetic neutrinos expected in the TD scenario are likely to exceed 10 14 GeV even at 
high redshifts of z > 100, collisions of these superhigh energy neutrinos with the cosmic background 
neutrinos should occur. The cosmic background neutrinos follow the Fermi distribution: 



(4) 
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Here k b b is the energy of the cosmic background neutrinos, T v is the present black body 
temperature (1.9 K), ks is the Boltzmann constant, and z is the redshift value at a given epoch. 

The collision probability between EHE neutrinos and the thermal background neutrinos is 
given by ( Yoshida 1994j ) 



dF , _ . 1 f , dcr , . f „. dn „ m , dkhh , , 

where is the collision angle between the EHE neutrino and the background neutrino, J5„ is the 
EHE neutrino energy at epoch with redshift z, s is the Lorentz invariant parameter which can be 
written in the Cosmic Ray Laboratory Frame (CRLF) as 

s = 2E u k bb (l + cos0), (6) 

a is the cross section for the interaction, and E rec is the recoil energy of collision particles such as 
electrons, muons, taus, quarks and neutrinos. Their decays or hadronizing process would produce 
neutrinos and further contribute to the neutrino cascade. The main interactions in the neutrino 
cascade are the following: a) s-channel Zq exchange, b)t-channel exchange, c) t-channel Zq 
exchange, d) interference between a) and b), e) interference between a) and c), and f) interference 
between c) and u-channel Zq exchange. The most important ones are channel a) which has 
resonance behavior, and channel b) which becomes important above the resonance energy for 
coupling of neutrinos with different flavors. The cross sections for these channels are respectively 
written as 

da ° 2s Ml r 5 l(i + cosr) 2 + 5 |(i-cosr) 2 j, (?) 

2' (8) 



dcos9* 4vr (s - Af 2 ) 2 + Mf r| 

da G 2 s Ml(l + cos6*) 2 



dcosfl* 4vr - cc .sfl*) + M 2 ] 

where G is the weak interaction coupling constant, 9* is the rotation angle of the collision in the 
center of momentum system (CMS), M z and M w are the masses of Zq and W^, respectively, 
r z is the decay width of Zq, and gi and gn are the left-handed and right-handed coefficients 
respectively. Though these two channels are the main contributors to the cascade, we take into 
account also the other channels c) through f). Since the energy of the recoil particles, E rec , is 
related to 9* through a Lorentz transformation, da/dE rec can be calculated from da/d9*. 

The number of produced neutrinos per unit length per unit energy is given by 

dNj^i f dn dF 

dE rec — — {E v> i, E rec ) — (E U) j,z), (9) 



dE u ^dL J dE v ^ ' dLdE,, 

where i and j denote the types of neutrinos like z/ e ,^ M , and v T , and dn/dE v j is the energy 
distribution of the secondary neutrinos produced by the recoil particle. For elastic collisions, 
these will be delta-functions, and for the produced muons and taus, they can be derived by 
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the decay matrix elements ( paisser 1990| ), Quarks fragment and produce jets of hadrons which 
emit neutrinos through their decays, so dn/dE V i can be gotten by convolution of the hadronic 
fragmentation spectrum with the hadron decay spectrum. The fragmentation spectrum can be 
written as ( [Hill et al. 19871 ) 

dN h 



/ _ 0.08 exp [2.6 y / ln(l/x)](l - xy[x\J\n(l/ x)}~ L . (10) 

Here, x = E/Ej e t is the energy of a hadron in the jet and Nh is the number of hadrons with 
fraction x of the energy in the jet Ej e t- We assume that 97% of the hadrons in the jet are pions 
( Bhattacharjee et al., 1992 ). 



Then we can write the transport equation of EHE neutrinos using Eq. (9) as follows ( Yoshida 



dN if „ ^_ Ni(E Vti ,z) , ^ [ E »,r w ATf ^ dN^iiE'^E^z) 



1994| ): 

+- fE r-[H (l + z)%E l/ji N i (E I/A ,z)]. i = v e ,u ll ,u T (11) 

Here Aj is the mean free path which is calculated by integrating Eq. (5) over E rec and s, Hq is the 
Hubble constant and the third term on the right hand side of Eq. (11) expresses the adiabatic loss 
due to the expansion of the Universe. Numerical solution of Eq. (11) gives the energy distribution 
of EHE neutrinos after their propagation. 

Figure |3| shows an example of the EHE neutrino energy distribution in the neutrino cascade. 
One can see the secondary bulk of lower energy neutrinos. This is mainly contributed by the 
hadronic jet which channel a) produces at around the Zq resonance energy (~ 3 x 10 15 GeV). You 
will see that this bulk enhances the final flux of neutrinos. 



4.2. Propagation of EHE neutrinos: Massive Neutrinos 

It is well known that invoking a small mass for cosmological background neutrinos might 
resolve the dark matter problem. Comparing the cosmological mass density and the neutrino 
density gives a constraint on the masses of stable neutrinos as ( |Bludman 1992 ). 

m v = 92n h 2 , [eV] (12) 

where £lo is the total cosmological mass density normalized by the critical density of the Universe, 
and the Hubble constant is Hq = lOO/i km sec -1 Mpc" 1 . The mass density is related to the 
Hubble constant Hq and the present age of the Universe io : 

{-q^T + (n -l)3/2 tan" 1 V^o - 1 Closed Universe, 
— q 1 1_ 1 — ( 1 _q°)3/2 tanh" 1 a/1 — Open Universe, (13) 
| Flat Universe 
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The age of the Universe can be estimated from measured ages of star clusters ( (Bandage and 



Cacciari 1990| ) to be t = 13 ~ 17 Gyr. Allowing 0.5 < h < 1.0, Eq. (13) and t give 

n h 2 < 0.038 ~ 0.255 (14) 
and the neutrino mass limitation is given by 

rn u < 3.5 ~ 23 eV. (15) 

Thus, for example, m Ve = ra Vli = m Vr = 1 eV are still allowed by this somewhat strict constraint. 
For simplicity in the following discussion, masses of neutrinos are assumed to be the same for all 
flavors. 

If the rest mass energy of cosmological background neutrinos is much higher than their black 
body temperature, their energy distribution is effectively monochromatic. Thus, the probability 
of collision between EHE neutrinos and the cosmic background neutrinos, Eq.(5), becomes 

dF , „ . . . o da , . 

■(E v ,z) ~ {1+zfno—- \ s= 2m„E„, (16) 



dtdE rec dE r 

where no is the number density of the background neutrinos in the present Universe (~ 108cm -3 
for v + v in each flavor). We get the mean free path by integrating Eq. (16): 

\{E U , z) ~ (1 + z)- 3 [noals^^EvV 1 ■ (17) 

Figure |] shows the mean free path obtained by Eq. (17). The Zq resonance effect on the cross 
section for channel a) (see Eq. (7)) creates the bump at around 

E v \resonance = M 2 j2m v = 4 X H) 12 ^)" 1 . [GeV] (18) 

lev 

Above this energy region, t — exchange (see Eq. (8)) become the dominant channels. 

Numerical solution of the transport equation, Eq. (11) , together with Eqs. (9), (16) and 
(17) give the EHE neutrino energy distribution propagating in the massive background neutrino 
field. In figure |5| is shown the energy distribution of EHE neutrinos propagating in the background 
neutrino field with m u = leV. The transition of the secondary neutrino energy distribution 
is found at around 10 13 GeV because the main contribution for the bulk of the lower energy 
neutrinos is the decay of pions in the hadronic jets which are emitted through the interaction of 
s — Zq exchange, while the secondary neutrinos with higher energies come mainly from the vv 
coupling through the exchange. 

Figure || shows the expected spectral shape from a single source of EHE neutrinos. The sharp 
dip appears because of the Zq resonance behavior as described above. Creation of this dip has 
been pointed out by Roulet ( [Roulet 1993 ) using simple analysis considering only the absorption 



effect by the background neutrinos. We also find a slight enhancement of the flux at lower energy. 
This enhancement created by the bulk of the secondary neutrinos will play an important role in 
studying the final flux of neutrinos from TDs. 
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4.3. Flux of EHE neutrinos 

The release rate of the X-particles by collapse of cosmic string loops or monopolonium 
( Bhattacharjee and Sigl 1995| ) is given by 



(19, 

where k is determined mainly from the abundance of the string loops in the lowest oscillation 
state or the abundance of monopolonium, and the X-particle mass, all of which are unknown. We 
get the final flux of EHE neutrinos after propagating in the cosmic background neutrinos using 
Eq. (1): 

w-r*%{m)'L dEkGl *--*'- t;),i *- ) - (20> 

The "green function" G is calculated by the transport equation, Eq. (11). Because the primary 
neutrinos are produced by decay of the pions in the hadronic jets from X-particle decay, the 
primary neutrino spectrum f(Et e ) is given by the convolution of the hadronic fragmentation 
spectrum, Eq.(10), and the pion-muon decay spectrum flYoshida 1994 , Bhattacharjee et al., 



1992] phattacharjee and Sigll995| ). We assume the jet energy is mx/2 and 97% of hadrons in 



the jet are pions as in the previous calculations (e.g., [Bhattacharjee et al., 1992] ). It should be 
noted that the EHE neutrino spectrum is not sensitive to the rate of annihilations in the early 
radiation-dominated universe [z > 10 4 ) since any neutrinos produced then would have cascaded 
to lower energies. It is sensitive, however, to the evolution described by equation (19) for z > 100 
both because of the high TD annihilation rate at early times and because the target background 
neutrinos then had higher energy and greater density. 

Because the value of k in Eq.(19) is uncertain, the observed flux of EHE cosmic rays 
sets the normalization. If the extragalactic component of EHE cosmic rays dominates above 
around 10 18 ' 5 ~ 10 19 eV where the change of the spectral slope has been observed ( |Bird et~aL 



1994, Nagano et al. 1992, Yoshida et al., 1995), it is consistent to assume that the cosmic rays 
above about 5 x 10 19 eV are mainly from the collapse of TDs. This assumption gives an intensity 
of protons from TDs which is lower than the observed flux over the whole energy range, and the 



superhigh energetic events observed by Fly's Eye ( Bird et al. 1995] ) and AGASA flHayashida et 



al. 1994) well beyond 10 20 eV can be explained as gamma rays from TD annihilations (Sigl et al. 



1994] |Bhattacharjee and Sigl 1995] |Chi et al. 1993|) . Here we use the normalization obtained from 
the assumption that the cosmic ray intensity at 5 x 10 19 eV observed by Fly's Eye ( Bird et al. 



1994|) is contributed by protons (and neutrons) from TDs. This normalization is constrained by 
their contribution to the diffuse gamma ray background for EL ~ 200 MeV, but not completely 
ruled out yet ( Sigl et al. 1995| ). The recent reevaluation of the gamma ray flux ( [Sigl et aL 



1996) demonstrates that the normalization is allowed for extragalactic magnetic field strengths of 



< 10 11 G, which are consistent with current estimates ( Kronberg et al. 1994| ). 



Figures and || show the expected flux from cosmic string loops or monopolonium for the 
case of massless neutrinos and those with leV mass, respectively. It is found that the secondary 
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neutrinos produced in the neutrino cascade enhance the intensity below ~ 10 20 eV. For the 
massive neutrino case, the sharp dip structure appearing in the spectrum from each source as 
shown in figure |6] suppresses the intensity at around the Z resonance energy expressed by Eq. 
(18). The flux at the lower energies, however, is enhanced by the cascading effect. 

The fluxes obtained here are much higher than the Greisen neutrinos discussed in section 2. 
The integral flux above 10 19 eV is listed in table 1. We should note that the flux estimate has 
an uncertainty of about a factor of 2 due to the poor cosmic ray statistics at 5 x 10 19 eV and 
possible systematics in the cosmic ray energy estimation ([Bird et al. 1994 |Yoshida et al., 1995| ). 



The flux depends on the assumption of the hadronic fragmentation spectrum because the hadronic 
jets play an important role in determining both the primary spectrum and the flux enhancement 



by neutrino cascading. It has been pointed out ( Chi et al. 1993 ) that the fragmentation spectrum 
could be steeper than Eq. (10), giving ~ E 1-1 ' 3 , and the final flux could be higher by more than a 
factor of 10. Thus the fluxes we present here could still be conservative. 

There is another potential TD which could be a source of EHE cosmic rays: the saturated 
superconducting cosmic string loops (SCS) ([Hill et al. 19871 ) . Although this hypoth esis might have 



several difficulties in producing EHE particles due to constraints by the low energy gamma ray 
backgrounds and the light element abundances in the Universe ( Sigl et al. 1995| ), strong evolution 



could give a high flux of neutrinos. The release rate of the X-particle is ~ 

t" 4 (~ t~ 3 for the 

ordinary string loops and monopoles as expressed in Eq. (19)), and the neutrinos at earlier epochs 
mainly contribute to the bulk of EHE neutrinos. A massive neutrino background would provide an 
especially interesting possibility. Figure ^ shows the expected EHE neutrino spectrum from SCS 
for a neutrino mass of 1 eV. The dip structure as seen in figure |6] is not smeared even in summing 
over the EHE neutrinos from SCS at different epochs according to Eq. (20), because of the strong 
evolution. The neutrino spectrum expected from this scenario would have this indicator of the 
neutrino mass. 



5. Detection of the EHE Neutrinos 

The fluxes calculated above are in the range 10 ~ 30vs km~ 2 yr~ 1 sr~ 1 above 10 19 eV and 
their detection requires a detector of huge aperture considering the low neutrino cross sections. 
The standard type underground neutrino telescopes, which observe long-range upward-going 
muons with effective areas of 0.1 km? or smaller, would not be capable of EeV neutrino detection 



( Gandhi et al. 1995 ). An alternative detection method is to search for extensive air showers (EAS) 
initiated by electrons produced by neutrinos through the charged current process u e + N — > e + X. 
Showers developing deep in the atmosphere must be produced by penetrating particles. The Fly's 
Eye experiment searched for such an event to get an upper bound on the EHE neutrino flux 



( Baltrusaitis et al., 1985 , Emerson 1992 ). The predicted fluxes are still lower than this bound by 
about a factor of 10. There are large new detectors for measuring EHE cosmic ray air showers 
using air fluorescence which are now under construction or development. These will have large 
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enough apertures to have the potential to search for EHE neutrinos (The HiRes collaboration 



19931 , [Teshima et al. 1992j ). These detectors can reconstruct the EAS development function 



of atmospheric depth with better than 30 g/cm 2 resolution and easily distinguish normally 
developing showers from deeply penetrating showers. 

The interaction length of cosmic ray hadrons and gamma rays is 50 ~ 100 g/cm? near 
10 19 eV. Thus the probability of these particles initiating EAS at deeper than 2000 g/cm 2 is less 
than 2 x 10~ 9 , so any shower starting that deep in the atmosphere would be a candidate neutrino 
event. The event rate of deeply penetrating showers (DPS) induced by EHE neutrinos is given by 

oo a u e 

^(>E Ue )=N A J dE' Ue J(E'j J dE e ^(E'j J dn(X n - 2000[g/cm 2 ])A(n\ (21) 

E Ue 

where a Ue is the cross section of the charged current process with the nucleon, E e is the energy of 
the produced electrons, Na is Avogadro's number, Xq is the slant depth of the atmosphere for the 
solid angle O, and A(Q) is the acceptance of the air fluorescence detector for deeply penetrating 
showers. 

The charged current cross section in this energy region has not been measured and could 
be uncertain. The observed small- x behavior of the QCD structure functions, however, gives an 
estimate of the cross section without large ambiguity and it turns out that the cross section is 



not saturated because of the evolution of the QCD structure functions (e.g., Mackey and Ralston 



1986] |Quigg et al.~1986| , Paisser and Stanev 198$ . Fi gure rL0j shows the calculated total cross 



section of v + N / ± + X. The increase in the cross section is favorable for EHE neutrino 
detection. 

The DPS acceptance, A(Q) in Eq. (21), is estimated by Monte Carlo detector simulations. 
For discussing the sensitivity for neutrino detection, it is useful to define the "effective aperture" 
which is folded with the cross section: 

D v {E Ve ) = N A J d E e ^(E Ue )T(E e ), (22) 
o 

T(E e ) = J dQ(X n - 2000[g / cm 2 }) A(Q, E e ). (23) 

T(E e ) represents the effective aperture column density (km 2 sr g/cm 2 ) of the detector for the 
DPS as a function of energy of the neutrino induced electrons. Then the event rate is estimated 
from Eq. (21) as 

oo 

dN 



dt 



(> E Ue ) = J dE' Va J{E Ve )D v {E Ve ). (24) 



Monte Carlo detector simulation can calculate the effective target volume T(E e ), but we shall 
give a rough estimate with a simple analytic calculation first. A phototube in the telescope is 
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useful in an event only if it collects more air fluorescence light emitted from the shower track than 
the fluctuation of night sky background light during its integration time t ga t e - The expected air 
fluorescence signal is given by 

N p h = Am ™~ N ^ ex.p(-r p /r )e eff r p AO (25) 

where r p is the shower's impact parameter, A m i r is the area of a mirror in the telescope, N e is 
the number of electrons in the shower cascade viewed by that phototube, ro is the extinction 
length of light due to the atmospheric scattering, e e ff is the fluorescence light yield from an 
electron (photons per meter), AO is the phototube pixel size and Q is the quantum efficiency of 
the phototube. The background light is given by 

N BG = n NBt gate A mir Q (AO) 2 (26) 

where n^B is the night sky photon intensity and t ga te is the gate time for collecting signal. Then 
the signal to noise ratio nth gives the threshold shower electron size for triggering a channel as a 
function of r p as follows: 



N e , th = n th 4nr p ex P (^Km/^% (27) 



1 V -^mtr^x 

To an accuracy of 35 %, this equation can be written as 

ro 



log 



n ( r \ ( e eff \ (R m ir\ I ( n NB \ ( t gate 



\8km J \Arn~ 1 ) \ lm J V \W 6 m~ 2 sr~ 1 fis~ 1 J \ 5fis 
where R m i r is the radius of the telescope mirror and Q is assumed to be 30 %. 



(28) 



The atmospheric slant width during which the shower cascade contains more electrons than 
this threshold size N e:t h can be considered as a target depth for showers to trigger a detector 
with almost 100 % efficiency. The longitudinal development of an electromagnetic air shower is 
described by the Greisen formula ( |Grciscn" T956| ) and we obtain the following expression for the 
target depth numerically: 

X} 00% =X t (N e >N etth ) = W0(-rj 2 -8r ] + 2) [g/cm 2 ] (29) 
rj = log(N etth ) - log (yj^) 

'100% : c 



Using Eq.(28), v can be written as a function of r p and thus X t is a function of E e and r p . 

r 1( 



We should remark that requiring X} 00% > leads to a maximum shower distance at which the 



telescope will trigger: 

4 



8km 



r rnax = ^ j _^0_ , - f ^ (30) 
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where 



log 



n t h 



( r o 



\8km ) \4m 



*eff 




n NB 



gate 



lm J V \ W b m 2 sr L us 1 J \ 5ns 

29km for the detectors under development 



(31) 



At W w eV with n th = 2 (2 a diviation) r™ ax 
operating in a desert atmosphere. 

Now we can calculate the effective aperture column density, T(E e ), as 
T(E e ) = 2tt [dcos9 [ dr p r p X™ 0% (r p ,E e ) [ d X . 



(32) 



Here we assume that our detector has azimuthally symmetric sensitivity. For integration over 
zenith angle 9, and azimuthal angle around the r p axis Xi we need to take into account the limit 
that the shower must be deeper than 2000g/cm 2 : 



Xq exp 



r p sin 9 sin x 
h 



> 2000 cos 9. 



[g/cm 2 



(33) 



Here Xq is the atmospheric vertical depth at the detector (~ 1000g/cm 2 ), and h is the atmospheric 
scale height (~ 7.5 km). After introducing some approximations, we finally obtained 



T(E e 



1.1 x 10° 



h 



8Akm 



8km 



f 2 -7x W 4 



h 



Akm 



[km 2 sr g/cm 2 



(34) 



Eqs.(22), (31), and (34) give the effective neutrino aperture for a detectors with mirror radius 
Rmir, and integration time t ga t e - 

We can confirm these estimates using the Monte Carlo program for the High Resolution Fly's 
Eye (HiRes), now being constructed at the Dugway, Utah ( The HiRes collaboration 1993j ). HiRes 
has two eyes separated by 12.6 km and 1° pixel size. We use the Monte Carlo code developed for 
estimate of the HiRes aperture taking into account the detector performance and atmospheric 
scattering ( pai 1993| ). 



Figure [Tl] shows the effective aperture for the HiRes detector obtained from both the Monte 
Carlo and our analytical form with R m ir = ^m, h = 7.5km, t ga t e = 5ns, and nth = 2. One finds 
that the Monte Carlo and the analytic methods are in good agreement. The neutrino aperture 
expands with energy both because of the increasing effective aperture column density for shower 
detection and because of the increasing interaction cross section, 

Beyond HiRes, the Japanese group is planning a large array of air fluorescence and Cherenkov 
telescopes named the Telescope Array ( [Teshima et al. 1992 ). We estimate the expected aperture 
for this experiment with possible parameters of R m i r = 1.5m and t ga t e = 350ns. It is also 



shown in figure 11. Though the final configuration and specification for the Telescope Array have 
not been determined yet and the aperture is somewhat uncertain, we should point out that an 
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air fluorescence detector could basically provide a neutrino detection sensitivity D u as high as 
~ 10 4 m 2 sr at 10 19 ey, which is determined mainly by the extinction length rp, the charged current 
cross section, and the solid angle for detection of the DPS rather than specific detector parameters 
like R mir and t gate - 

The event rate for neutrino induced showers is estimated by Eq.(24). Table 2 gives the 
expected event rate in HiRes and the Telescope Array for the possible EHE neutrino flux of 
the TD models. Figure ^ shows the detection probability during 10 years of running. It is 
found that the SCS model gives a measurable flux of neutrinos. For the ordinary cosmic string 
or monopole scenario, whether or not EHE neutrinos are detectable depends on the X-particle 
mass, which determines the enhancement factor of the fluxes due to neutrino cascading. If the 
mass of X-particles is heavier than 10 16 GeV, we have a chance to detect EHE neutrinos with air 
fluorescence detectors during 10 years of observation. 

An uncertainty in this estimate arises from uncertainty of the primary neutrino flux estimate 
(at least a factor of 2), and from uncertainty in estimating cross sections for the interactions 
of EHE neutrinos with nucleons. Gandhi et al recently reevaluated the cross sections based on 
improved knowledge of parton distributions measured by the electron-proton collider HERA 
(Gandhi ct al. 1995). The dashed curve in figure [H] shows their new estimate. It is larger than the 
calculated value based on the previous estimation (e.g., Mackey and Ralston 1986)) by 20 % at 1 
EeV, 40 % at 10 EeV, and 70 % at lOOEeV. This reevaluation, however, would not significantly 
affect the event rates listed in table 2 since statistical fluctuations and other uncertainties (for 
example, the fragmentation spectrum equation (10)) would dominate an uncertainty due to 
ambiguity in the cross section. 



6. Backgrounds 

Neutrino interactions in the atmosphere are expected to produce at most a few detections of 
deeply penetrating air showers over the course of a long experiment. It is therefore important to 
study whether other phenomena could produce similar deep showers. One possibility is a DPS 
resulting from the decay of an EHE r lepton deep in the atmosphere. Another possibility is a DPS 
due to an EHE 7-ray bremsstrahlung by a muon. In either case, the DPS would be associated 
with a higher energy primary shower which produced the r or [i. Although the primary shower 
may be far away, it should still be evident as an EHE shower either by its fluorescence light or its 
Cherenkov light. The danger of mistaking the secondary shower for a neutrino event arises only 
if distant clouds block all evidence of the primary shower without obscuring the DPS secondary 
shower. 

To be conservative, we will here assume that the primary shower is not detected when a deep 
EHE shower is produced by tau decay or mu bremsstrahlung. We evaluate an approximate upper 
limit for the rate of detectable deep secondary showers of that type caused by cosmic ray air 
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showers in the atmosphere. Several factors combine to make this rate small. First of all, air shower 
production of EHE taus and muons is small enough that their intensities are low compared to 
the intensity of equal-energy primary cosmic rays. Secondly, they can only occur at zenith angles 
between 70° and 90°. Otherwise there is not enough slant depth along the axis for the secondary 
shower to start deeper than 2000 g/cm 2 and develop above ground level (taken to be at 860 g/cm 2 
vertical depth for the analysis here). Thirdly, on these long slant shower axes, the probability is 
small that the secondary shower will develop within the detectable volume of atmosphere. The 
secondary shower is likely to start too soon or too late. Taken together, we find that these factors 
cause the detection rate for deeply penetrating secondary showers to be almost negligible for the 
detectors now being built or planned. EHE neutrino astronomy with air fluorescence detectors is 
limited by the signal, not by the background. 

Extremely high energy taus and muons may be produced in a shower either from the weak 
decay of a heavy quark, from the decay of a (real or virtual) W or Z, or from Drell-Yan processes. 
Being weak processes the contributions from the weak boson and Drell-Yan production are 
negligible compared to that from the production of heavy quarks via the strong interaction. The 
production rate of EHE taus and muons have therefore been estimated from tic bb cross-sections, 
using a first-order perturbative QCD calculation as implemented in the PYTHIA event generator 
( Sjostrand 1992|) , and using the default fragmentation parameterization. Each produced bottom 
or charm particle is allowed to decay (rather than interact), and the numbers of resulting taus 
and muons are counted at each energy. In a realistic simulation, the hadronic cascades should be 
simulated along sampled trajectories through the atmosphere, and many of the bottom and charm 
particles would interact instead of decaying. That would degrade the bottom and charm to lower 
energies, and the effect would be to reduce the overall intensity of taus and muons at all energies. 
By allowing each produced bottom or charm particle to decay, we are overestimating the tau and 
muon intensities. The resulting estimate for deep secondary showers is therefore only useful as an 
upper limit. The true rates should be lower. 

These (upper limit) intensities of taus and muons depend on the cosmic ray spectrum. If the 
spectrum is a power law, then the rates depend on the spectral index and the maximum cosmic ray 
energy (spectrum cutoff). The intensity is larger for flatter spectra and for larger cutoffs. Table 
3 displays results for some different assumptions about the cosmic ray differential spectral index 
and cutoff energy. For each case, the integral intensity of cosmic rays above 10 19 eV is normalized 
to the observed value of 0.5/(km 2 sr yr). It should be emphasized that the spectrum hypotheses 
are for the cosmic ray spectrum at Earth, not the source spectrum. The spectral hypotheses in the 
table are probably unreasonably hard, in view of the pion photoproduction energy losses. The r 
and fi intensity upper limits would be lower for more plausible softer spectra. 

For the case of a secondary tau lepton, it is also necessary to consider the probability that it 
will decay deep enough to be confused with a neutrino shower (>2000 g/cm 2 ) but early enough 
for its shower development to be above ground and within the fiducial volume which is consistent 
with equation (30). This probability is small for the high energy taus. The mean decay length 
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is 50 km at 1 EeV, 500 km at 10 EeV, etc. The probability depends on the trajectory and the 
r point of production as well as the r energy. We have numerically integrated over all viable 
trajectories, assuming the r is produced at 500 g/cm 2 slant depth. (This is deeper than expected 
for production of a high energy r. It makes the detection volume closer than expected to the tau 
production point, consistent with our pursuit of an upper limit.) We also assume that all of the 
tau's energy goes into the secondary air shower. In reality some tau decays produce only muons 
and neutrinos, yielding no air shower. In the other decays, at least some of the energy is taken by 
neutrinos, so we are also overestimating the detectable deep tau shower rate by not simulating the 
tau decays in detail. Results for various cosmic ray spectra are given in table 3. The tabulated 
upper limit for detected deep r showers is for a 10-year experiment with a 10% duty cycle. Even 
for the hardest spectrum, the expected number of deep showers due to taus is less than 0.1. 

The deep muon rates in table 3 result from numerically integrating over possible trajectories, 
as in the deep r rate calculation. In the muon case, the computation includes the probability that 
the muon will produce some bremsstrahlung gamma ray whose shower starts deeper than 2000 
g/cm 2 and is detectable in the fiducial volume. The expected number of deep showers due to 
muon bremsstrahlung is even less than what is expected from tau decays. 

If a fluorescence detector experiment were expected to detect many deep secondary particle 
EHE air showers, then there would be concern that one or more of them might be mistaken for 
a neutrino event due to the primary shower being obscured by clouds. The estimates here are 
unrealistically high and serve only as upper limits. They show that no detections of secondary 
deep showers are expected. With or without clouds, therefore, there should be no such background 
to interfere with the identification of neutrino events. 

The search for deep secondary showers is of interest in its own right. The estimates here 
show that none are expected. The detection of deep secondary showers would be evidence for 
an enhancement of bottom or charm production. It is therefore pertinent to search for deeply 
penetrating EHE subshowers which occur in association with distant primary air showers of higher 
energy. 

7. Conclusions 

We have analyzed the production of EHE neutrinos and estimated their flux based on several 
different models. The most certain process for producing such neutrinos - the decay of pions 
photoproduced by EHE cosmic rays propagating in the microwave background radiation - is not 
able to give rise to measurable EHE neutrino fluxes unless sources with extremely strong evolution 
emit EHE cosmic rays with very hard spectra. In that case the neutrino flux would exceed the 
cosmic ray flux by a factor of more than 100 at around 10 19 eV. However, it is hard to reconcile 
such strong evolution with deep sky x-ray and radio surveys. 

An allowed mechanism which would result in copious EHE neutrinos is the annihilation 
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or collapse of topological defects (TDs) such as monopoles and cosmic strings, which could 
also account for the highest energy cosmic ray events recently detected well beyond 10 20 eV. 
The energies of EHE neutrinos produced by this process range up to the GUT energy scale 
(~ 10 16 GeV) and their interactions with cosmic background neutrinos initiate significant neutrino 
cascading. The cascade neutrinos enhance the flux at earth, and we have derived the expected 
flux, including this enhancement factor, both for the assumption of stable massive neutrinos and 
for the assumption of massless neutrinos. We find that the neutrino intensities would dominate 
over the observed cosmic ray intensity at 10 19 eV by a factor of 30 ~ 100 with reasonable 
assumptions concerning the mass of the X-particles and neutrinos and using the expected evolution 
of monopoles and cosmic strings. Heavier X-particles would produce a higher neutrino flux due to 
increased enhancement from neutrino cascading. If neutrinos have limited mass, the Zq resonance 
effect on the neutrino cross section results in a slight suppression of the neutrino intensity at the 
resonance energy together with enhancement of the fluxes at lower energies. For TDs with stronger 
evolution such as the saturated superconducting cosmic string loops, the resonance effect would 
create a clear dip structure in the neutrino spectrum after propagating in a massive neutrino 
background. Therefore the neutrino spectrum expected from the TD scenario would contain an 
indicator of the neutrino mass. 

The large new detectors being built to observe EHE cosmic ray air showers by air fluorescence 
have the potential to detect EHE neutrinos as anomalous showers which start deep in the 
atmosphere. We have evaluated the sensitivity of these detectors for discovering an EHE neutrino 
flux. The conclusion is that one of the new cosmic ray detectors could detect a few EHE neutrino 
showers during 10 years of observation if the topological defect scenario is correct. A potential 
background of deeply penetrating showers could arise from secondary showers which result from 
r lepton decays deep in the atmosphere or from EHE gamma ray bremsstrahlung by muons. In 
either of these cases, it should be possible to recognize the higher energy primary shower which 
gave rise to the tau or muon. Such secondary showers need not be a genuine background for 
neutrino shower detection. Even if the primary showers were not detectable, our estimated upper 
limits for DPS event rates produced by high energy secondary rs and fis are found to be more 
than two orders of magnitude lower than the neutrino event rates expected from the TD scenario. 
Therefore EHE neutrino astronomy with air fluorescence detectors is limited by the signal, not by 
the background. A search for EHE neutrinos is indeed a meaningful test of the TD hypothesis. 
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Fig. 1. — Energy spectrum of the Greisen neutrinos under several assumptions of m and z, 
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Fig. 2.- 



Energy spectrum of the Greisen neutrinos for the strong source evolution case. 
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Fig. 3. — The energy distribution of EHE neutrinos in the cascade after propagation through IGpc. 
Primary input spectrum is monochromatic energy distribution of 10 16 GeV of muon neutrinos. 
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Fig. 4. — The mean free path of EHE neutrinos in the massive cosmological neutrino backgrounds 
at present epoch. 



10 12 14 16 10 12 14 16 

LOG(Energy) [GeV] LOG(Energy) [GeV] 

Fig. 5. — The energy distribution of EHE neutrinos after propagation through lGpc in the 
background neutrinos with m v = 1 eV . Primary input spectra are monochromatic energy 
distribution of 10 15 GeV and 3 x 10 12 GeV of muon neutrinos respectively 




Fig. 6. — The energy spectra of EHE neutrinos propagating in the massive cosmological background 
neutrino field. The primary spectra are assumed to be ~ E~ 2 . 
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Fig. 7. — The differential flux of the EHE neutrinos from the collapse of the cosmic string loops or 
monopolonium. All neutrinos are assumed to be massless particles. The dashed curves show the 
results calculated if only the energy loss due to the expansion of the Universe is considered. 
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ENERGY [GeV] 

Fig. 8. — The differential flux of the EHE neutrinos from the collapse of the cosmic string loops 
or monopolonium. Masses of all neutrinos are assumed to be 1 eV. The dashed curves show the 
results calculated if only the energy loss due to the expansion of the Universe is considered. 
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Fig. 9. — The differential flux of the EHE neutrinos from the saturated superconducting cosmic 
string loops. Masses of all neutrinos are assumed to be 1 eV. 
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Fig. 10. — The cross section of the charged current process v + N — ► + X calculated by Mackey 
and Ralston. The dashed curve was obtained by the recent reevaluation of parton distribution 
measured by HERA. 
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Fig. 11. — The effective aperture for the EHE electron neutrinos expected for the HiRes and the 
Telescope Array detector performances. The real curve shows the results obtained by the HiRes 
Monte Carlo simulation, the dashed curve was obtained by the analytical formula with the HiRes 
detector parameters, and the dotted curve corresponds to the aperture of the Telescope Array 
which was also calculated by the analytical formula. 



-32- 
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Fig. 12. — Detection probability of the EHE neutrinos during 10 years observation of the HiRes 
(the real curves) and the Telescope Array (the dotted curves). 10 % duty cycle is taken into 
account. The plotted detection probability is the Poisson probability of detecting at least one 
neutrino shower, given the calculated expected number. 
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m x = 10 16 GeV m x = 10 15 GeV m x = 10 15 GeV 
m v — eV m v — eV m v = 1 eV 


Observed 
Cosmic Rays 


V V 


1.4 x 1CT 12 4.4 x 1CT 13 3.8 x 1CT 13 
2.8xl0~ 12 8.9 x 10~ 13 7.7xl0~ 13 


1.4 x 1(T 14 



Table 1: Integral Neutrino Flux above 10 19 eV for the different assumption [m 2 sec 1 sr 1 ]. 





m x = 10 16 GeV to* = 10 15 GeV m x = 10 15 GeV SCS 

= eV m v = eV m v — 1 eV m u = 1 eV 


E > 10 18 eV 
E > 10 19 eV 


0.38(HR) 0.79(TA) O.ll(HR) 0.23(TA) O.ll(HR) 0.22(TA) 5.93(HR) 15.6(TA) 
0.23(HR) 0.41(TA) 0.08(HR) 0.14(TA) 0.07(HR) 0.12(TA) 0.99(HR) 1.77(TA) 



Table 2: The expected event rate by the HiRes and the Telescope Array induced by the EHE 
neutrinos from the TDs collapse [/5yr]. 



Differential 


Energy 


It/Icr 


Ifi/ICR 


Detectable 


Detectable 


spectral index 


Cutoff 


(E > 1 EeV) 


(E > 1 EeV) 


deep t showers 


deep fi showers 




(EeV) 


(upper limit) 


(upper limit) 


(upper limit) 


(upper limit) 


3.0 


1000 


7.5 x 10~ 8 


2.9 x 10~ 7 


2.9 x 10~ 4 


3.5 x 10~ 7 


2.0 


1000 


5.3 x 10~ 6 


2.3 x 10~ 5 


2.2 x 10~ 3 


9.7 x 10~ 6 


2.0 


10,000 


8.2 x 10~ 6 


3.7 x 10" 5 


3.4 x 10~ 3 


2.6 x 10" 5 


1.3 


1000 


1.6 x 10~ 4 


7.2 x 10~ 4 


1.6 x 10~ 2 


1.1 x 10" 4 


1.3 


10,000 


7.4 x 10~ 4 


3.5 x 10" 3 


6.6 x 10~ 2 


8.1 x 10" 4 



Table 3: Upper limits for the expected number of deep secondary showers during a 10-year 
experiment with 10% duty cycle (so 1 year of on-time). The first two columns describe the 
hypothetical primary cosmic ray spectra at Earth. The third column is the (upper limit) integral 
intensity of produced taus above 1 EeV, normalized by the integral intensity of primary cosmic 
rays. The fourth column is the corresponding quantity for muons. The last two columns give upper 
limits for the numbers of expected deep EHE secondary showers detected as the result of tau decays 
and muon bremsstrahlungs, respectively. 



